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In a first approximation, known as tiie adiabatic process, the direction of the spin poiarization of 
currents is parallel to the local magnetization vector in a domain wall. Thus the spatial variation of 
the direction of the spin current inside the domain wall results in an adiabatical spin transfer torque 
on the magnetization. We show that domain wall motion driven by this spin torque has many unique 
features that do not exist in the conventional wall motion driven by a magnetic field. By analytically 
and numerically solving the Landau-Lifshitz-Gilbert equation along with the adiabatic spin torque 
in magnetic nanowires, we find the domain wall has its maximum velocity at the initial application 
of the current but the velocity decreases to zero as the domain wall begins to deform during its 
motion. We have computed domain wall displacement and domain wall deformation of nanowires, 
and concluded that the spin torque based on the adiabatic propagation of the spin current in the 
domain wall is unable to maintain wall movement. We also introduce a novel concept of domain 
wall inductance to characterize the capacity of the spin-torque induced magnetic energy stored in a 
domain wall. In the presence of domain wall pinning centers, we construct a phase diagram for the 
domain wall depinning by the combined action of the magnetic field and the spin current. 

PACS numbers: 75.75.+a, 75.30.Ds, 75.60.Ch, 72.25.Ba 



I. INTRODUCTION 

The subject of current-induced magnetization reversal 
(CIMR) in magnetic multilayers has received consider- 
ably interest recently H S S 0, 0, . From fundamen- 
tal point of view, this topic introduces a new interaction 
between non-equilibrium conduction electrons and local 
moments, and the physics of this new interaction has not 
been well explored. From the application point of view, 
the current-induced magnetization reversal opens a novel 
way to control and manipulate the magnetization dynam- 
ics that is one of the central issues in modern magnetic 
technologies. In conventional magnetic devices, the mag- 
netization direction is controlled by external magnetic 
fields generated by a current or by a permanent magnet. 
These magnetic fields are usually spread into a relatively 
large spatial extent. The CIMR can be confined to an 
exact spatial region where the current flows; this prop- 
erty is very attractive for magnetic nano-devices, e.g., 
magnetic random access memory. 

The physics of the current-induced magnetization re- 
versal involves interplay between non-equilibrium con- 
duction electrons and local magnetization. Namely, the 
spin angular momentum carried by spin polarized elec- 
trons is transferred to the local moment through the ex- 
change interaction. This is somewhat analogue to the 
well-known phenomenon of electromigration 7] where 
(impurity) atoms move along the direction opposite to 
the direction of the current. The origin of the electromi- 
gration is mainly from the transfer of linear momentum 
of non-equilibrium conduction electrons to atoms, result- 
ing a "wind force" which drives migration of the atom. 
For magnetic materials where the electric current is spin 
polarized, it is possible that the spin-polarization of non- 
equilibrium conduction electrons changes its orientation 
along the direction of the current, and thus the local- 




FIG. 1: Cartesian and polar coordinate systems 

ized magnetic moment receives a spin torque if the con- 
duction electron continuously transfers their spin angular 
momentum to the local moment in the steady state of the 
current flow. This non-equilibrium conduction electron 
induced spin torque can alter magnetization structures, 
drive magnetization dynamics and even create the do- 
main wall motion. Up until now, the spin transfer model 
has been mainly on a trilayer structure where two fer- 
romagnetic layers are separated by a spacer layer. Typ- 
ically, one chooses the spacer layer thick enough to re- 
duce the magnetic coupling between the two ferromag- 
netic layers and thin enough to minimize the spin-flip 
scattering in the spacer so that the spin transfer torque is 
maximized. Several theoretical models 
have been put forward to formulate the spin torque in 
terms of geometric and materials parameters. While it is 
still debatable on the correct microscopic picture of the 
spin torque, all theories contain a spin torque in the form 
of 

= -^Mx (MxMp) (1) 

Ms 

where M is the magnetization of the free layer, Mp is the 
unit vector along the direction of the magnetization of the 
pinned layer, is the saturation magnetization, and 
aj is a model-dependent parameter and is proportional 
to the current density. In general, the magnitude aj of 
the spin transfer torque is also a function of the angle 
between the pinned and free layer magnetization vectors. 
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Several theoretical attempts 0,0,0,113 were made 
in understanding the unique features displayed in the 
spin torque, Eq. (1). Here we just list three of them 
below: 1) the spin torque may stabilize steady preces- 
sional states, i.e., the magnetization does not converge 
to a meta-stable static state even at zero temperature; 
2) the spin torque may create a new stable configuration 
of the magnetization, e.g., an out-of-plane magnetization 
direction can be a new solution of meta-stable states; 3) 
the spin torque modifies the effective energy barrier or 
temperature in a significant and unique way. 

All of the above novel properties are derived based on 
the specific form of the spin torque, Eq. (1). Equation 
(1) has assumed that both free and pinned layers have 
uniform magnetization within the layers, at least along 
the direction of the current. In a typical ferromagnet, the 
magnetization is rarely uniform and the dynamical pro- 
cess of magnetization is not coherent rotation in general. 
It would be interesting to find what the spin torque does 
to magnetization dy namics of non-single domain ferro- 
magnets. Berger [l^l introduced the "domain drag force" 
by considering the spin torque in a single layer system 
where the magnetization is not uniform along the cur- 
rent. He argued, based on his intuitive physics picture, 
that the current can drag the domain wall along the path 
of the current flow. Bazaliy et al. proposed a spin torque 
in a ferromagnet within the ballistic transport model for 
half-metallic materials H^]. Most recently, we have gen- 
eralized the spin torque for any diffusive transport ferro- 
magnet 13, 



Tst = -bjM X [M X (je • V)M] 



(2) 



where M = M/M„ je = je/je and hj = Pje^iB/eM„ P 
is the spin polarization of the current, /^s is Bohr mag- 
netron, and je is the electric current density. The above 
torque has identical form as that of Bazaliy et al [T^. 
However, we derived the above formulation on a general 
property of ferromagnetic materials: in a ferromagnet, 
the spin polarization of the current is always along the 
direction of the local magnetization vector, i.e., the trans- 
verse component of the spin current can be neglected 
jlOl [TIL ITsj . Thus one may define a spin current tensor 
J = {fJ'B/eMs)Pie (8) M(r). The spin transfer torque 
is defined as the rate change of the angular momentum 
of conduction electrons that are transferred to the local 
magnetization, Tgt — dm/dt — dm/dt + V-J. In the 
steady state of the current fiow, dm/dt — and thus we 
find the above spin torque can be written as Eq. (2) if 
we use the fact that the magnitude of the local magneti- 
zation vector is a constant. 

While the magnitude of the spin torque in spin valves, 
Eq. (1), depends on many unknown parameters such as 
interface spin-dependent resistance, the magnitude of the 
adiabatic spin transfer torque, bj in Eq. (2), can be 
accurately estimated for many ferromagnetic materials. 
bj, which has a dimension of velocity, is determined by 
two materials parameters: and P; these parameters 
have already been determined experimentally. In table 



I, we list the values of bj for a selected set of mate- 
rials for a spin current density je = IQ'^ A/ cm? . The 
half-metals, Cr02 and Fe20^^ has a maximum spin po- 
larization (P=l) and has low saturation magnetization, 
and thus the spin torque is larger than that of transition 
metals. 

Table I The values of bj for some materials for 

je = lO^A/cTO^ 



Nanowire 


Ms{A/m) 


P 


bj(m/s) 


Fc 


17.18 X 10=' 


0.5 


1.35 


Co 


14.46'' 


0.35 


1.41 


Ni 


4.9 X 10-' 


0.23 


2.7 


Permalloy 


8 X 10'' 


0.7 


5.1 


7 - Fe20-j, 


4.14 X lO'' 


1.0 


14.0 


CrOa 


3.98 X lO'' 


1.0 


14.6 



The paper is organized as follows. In Sec. II, we estab- 
lish the equation of motion of a domain wall by using the 
generalized Landau-Lifshitz-Gilbert equation along with 
the adiabatic torque, Eq. (2). We write the equation for 
nanowires. In Sec III, we propose a trial function to an- 
alytically obtain solutions for domain wall dynamics. In 
particular, we address the spin torque induced domain 
propagation and distortion. Then a detailed comparison 
between numerical results and analytic solutions is made. 
Finally, we summarize our major findings in Sec. V. 



II. MODEL SYSTEM 

We assume the current fiows in x-direction, along the 
long length of a wire, i.e., je = e^. By placing this into 
Eq. (2) and by treating Eq. (2) as an additional torque 
on the standard LLG equation, we write the generalized 
LLG equation in the presence of the spin torque as. 



dm 

~dt 



-7M X H, 



^aM X 



dm 

~dr 



bjM X I M X 



dM 

dx 



(3) 

where 7 is the gyromagnetic ratio, and He// is the ef- 
fective magnetic field including the external field, the 
anisotropy field, magnetostatic field, and the exchange 
field, and a is the Gilbert damping parameter. 

We now specify the geometry of the wire and the effec- 
tive field entering Eq. (3). As shown in Fig. 1, the x-axis 
is taken as the easy axis as well as the direction of the 
external field Hf-xt- The width and the thickness of the 
wire are D and d, respectively. A current is along the x 
direction. We explicitly write the effective field, 

H Al 2yl 
He// = ^/ "' e^ + -TT^V^M - AirRFe, + He^te^ (4) 
Ms 

where Hk is the anisotropy field, A is the exchange con- 
stant, and AttMz is the de-magnetization field. 

To analytically gain insight on the domain wall dynam- 
ics, we first consider that the magnetization M varies 
only in the direction of x-axis. By placing Eq. (4) into 
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(3) , the LLG equation can be conveniently written in the 
polar coordinate, 



1- Q, sm = 7 2 cos 6*-^ — + sm 9 — ^ 



dt 



dt 



M, 



dx dx 



de 



-477rMs sin0sin(/?cos(/5 + 6j— — (5) 

ax 



89 
"di 



-^Hk sin 9 cos 9 — ijirMs sin 9 cos 9 sin^ ip 




--fHext sin 9 - bj sin 6*— (6) 

where 9 represents the angle between the magnetization 
vector and x-axis, and if is the out-of-plane angle of 
the magnetization vector projected in the yz plane, i.e., 
Mz — Ms sin 9 sin (p. 

Our goal is to solve the magnetization vector as a func- 
tion of position x and time t from the above equations 
of motion. In the absence of a current and/or an ex- 
ternal magnetic field, we consider that the domain wall, 
separated by two head-to-head domains along the wire 
direction, is a Neel wall whose magnetization stays in 
the plane of the wire, i.e., — 0, = MsC0s9, 
My — Ms sin 9 cos (p, where 



(^ = 0; 



= 2tan ^exp(x/Wo) 



(7) 



and Wq = \/2A/ HkMs is the domain wall width. At 
t = 0, an electric current /an external field is applied. We 
should determine the motion of the wall at i > from 
Eqs. (5) and (6) by utilizing the initial condition, Eq. (7). 



III. ANALYTICAL RESULTS 

The non-linear partial-differential equations, Eqs. (5) 
and (6), are difficult to solve. Here wc follow Walker's 
analysis of domain wall motion by introducing a trial 
fimction !2(M. 



LP = ip{t); 



In tan ■ 



cit) 



w(r)dr . (8) 



The first equation assumes that the projection of the 
magnetization vector in the domain wall on the yz plane 
is independent of the position. One should note that the 
spatial independence of (p{t) does not mean a uniform 
out-of-plane component because M^ — sin 9 sin ip and 9 
varies spatially in the domain wall from 9 — to n. The 
second equation in Eq. (8) postulates that the domain 
wall shape remains a standard Neel- wall form except that 
the width of the wall W{t) = c{t)^^ varies with time and 
the wall moves with an instantaneous velocity v{t). As 
we will show later that this form of the solution is indeed 
a correct solution as long as the spin torque bj and the 
external field Hg^t are small. 



By placing this trial function into Eqs. (5) and (6), and 
by utilizing the following identities. 



dip d'^ip 
'd^~~dx^ 







961 
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— — = c(t) sin6'; ^TT ^ '^(^) sin^cos 



dx 



dx"^ 



do , , . 
— = ci[x, t) smt 
dt 



and defining a function ci{x, t), 

C^{X, t)^'^[x~ v{T)dT ] - C{t)v{t) (9) 



we find 



dip 



ci{x, t) + a— = —A-K^Ms ninLpcosip + bjc(t) (10) 
dt 

. dip 
-aci[x, t) + — - ^Hext = 
dt 

7 (^Hk ~ J^c{tf + ^T^Ms sin^ i^ cos 9 (11) 

Strictly speaking, Eqs. (10) and (11) can not be valid 
because ci(x,t) defined in Eq. (9) is linearly proportional 
to X but there are no other terms in Eq. (10) and (11) 
containing spatial variable x. Thus, the trial function of 
the form of Eq. (8) is only possible when one discards the 
spatial dependence of ci. Equivalently, the first term in 
the right side of Eq. (9) must be much smaller compared 
to the second term so that ci « —c{t)v{t). While one 
would not determine the smallness of this first term a 
prior, we will later confirm that c{t) only shrinks by a 
few percent during the domain wall motion and it varies 
slowly with time so that dc{t)/dt is indeed small. Now 
the left-hand side of Eq. (11) becomes ^-independent and 
thus the prefactor in front of cos 9 at the right-hand side 
of Eq. (11) must be identically zero, 

2A 

Hk c{tf + AttMs sin^ ip = (12) 

Ms 



(13) 



and 

ac(t)v{t) + — = -yHext ■ 
dt 

By placing Eq. (13) into (10), wc have 



(1 + ~ "jHext + Oibjcit) — 4na^Ms sin ip cos ip . 

dt 

(14) 

Equations (12) and (14) are the ordinary first-order dif- 
ferential equations that determine the domain wall width 
c^^{t) and the rotation of the domain wall plane (p{t) 
subject to the initial values Lp{t = 0) = and c^^{t — 
0) = I/Wq. Once these two equations are solved, we can 
obtain the velocity of the domain wall from Eq. (13) 



vit) 



^{aHext + AttMs sin ip cos (p) 
(l + a2)c(i) 



1 + 



(15) 
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FIG. 2: The velocity v, the wall width W, the wall angle 
ifi and the displacement function of time t for differ- 

ent spin torques and for H^xt ~ OOe. The parameters are 
4ttMs = 1.8 X 10* Oe, Hk = 500Oe, Ms = 14.46 x lO^A/m, 
^ = 2.0 X 10"" J/m and a = 0.02. 



B. Domain wall distortion 

In our analytical theory, domain wall distortion is char- 
acterized by two parameters, ip(t) and c~^{t). The former 
describes the out-of-plane component of the magnetiza- 
tion and the latter is the time-dependent domain wall 
width. When t oo, the distortion is maximum. By 
placing dip/dt = 0, we find 1^3(00) from Eq. (14) in the 
absence of the external field. 



sin 2(p(oo) = 



bjc{co) 



2A 



2A 



sin^ (fiioo) 



fl7) 



2j7tMs 277rMs 

where the last identity is from Eq. (12). When bj 
(<C 7Wov'47rAfsiJif ) is small, sinip ip and sin^ « 0, 
Eq. (17) is simplified as. 



AttjMsWo 



(18) 



By inserting this into Eq. (12) and again by assuming bj 
is small, we find the maximum reduction of the domain 
wall width is 



1 



W{co) _ ^ 
Wq ~ ~ At: MsHkI^ 



(19) 



In Fig. 2 (b) and (c), we show the domain wall distortion 
during the application of the current. 



The numerical solutions of Eq. (12) and (14) are repre- 
sented in Fig. 2 where the domain wall velocity, displace- 
ment, and distortion as a function of time are shown. In 
the following subsections, we discuss these quantities. 



A. Domain wall velocity 

At t=0, the wall is a standard in-plane Neel wall, i.e., 
tf{t = 0) = 0. From Eq. (15), we have readily seen that 
the velocity is —bj/{l-\-a^) at the initial apphcation of 
the current if no external field is applied. On the other 
hand, for t — > 00, dip/dt = 0. We immediately find, from 
Eq. (13), 



u(oo) 



c(oo)a 



(16) 



this result is same as that of Walker |20j in the absence of 
the spin torque. We conclude that the terminal velocity 
is independent of the spin torque at all! This implies 
that the spin current alone, i.e., no applied field, can not 
move the domain wall to a large distance. In Fig. 2(a), 
the velocity of the domain wall during the application of 
the current is shown. It is noted that the domain wall 
motion stops at a fraction of a nanosecond. 



C. Domain wall displacement 

As the domain wall motion eventually stops even for 
the perfect wire in the absence of the external field, it 
is interesting to see how far the domain wall moves, i.e., 
what is the maximum distance a domain wall can travel 
before it stops? We estimate this displacement by inte- 
grating the velocity and utilizing Eqs. (12) and (13) 



vdt 



1^(00) 



dip 
ac 



(20) 



Notice that the displacement is inversely proportional to 
the damping constant. In Fig. 2(d), we show the dis- 
placement as a function of time. 



D. Inductance of domain wall 

We have seen that domain wall motion driven by the 
current is so much different from that driven by an ex- 
ternal field. In the latter case, the initial velocity upon 
the application of the external field is small and the ve- 
locity is gradually increasing until a saturation velocity 
is reached; after that, the domain wall moves with a con- 
stant velocity. We can understand this field-driving wall 
motion rather straightforwardly: the wall motion along 
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the direction of the field is to reduce the Zeeman energy; 
the rate of the Zeeman energy reduction equals to the 
energy damping in a constant moving domain wall. For 
the current- driving domain wall motion, the wall motion 
does not decrease the energy because the spin torque has 
no cifect on the; uniformly magnetized domains. Instead, 
the application of the spin torque introduces energy to 
the wall so that distortion of the wall occurs. The con- 
sequence of the wall distortion is to generate an energy 
damping mechanism so that the pumping of the energy 
by the spin torque can be compensated by the damp- 
ing. Once the dynamic equilibrium between pumping 
and damping is reached, the domain wall motion stops 
and the distortion is the maximum. To quantitatively 
see this interesting feature, we consider the wall energy 




a=0.02, H„=0Oe 



(b) 



E 
E 



300 600 900 300 600 900 
Spin torque bj(m/s) Spin torque bj (m/s) 

FIG. 3: Maximum wall width narrowing AW/Wo and the 
maximum angle as a function of the spin torque. 



E 



2TrMt 



K 



2M., 



Mi 



A 



dM 



dx 



dx (21) 



where the first term is magnetostatic energy, the second is 
the anisotropy energy and the final term is the exchange 
energy. The rate of the energy change can be obtained 
via dE = —He// • dM. and use the LLG equation. 



rH 



eff 



dE 

dt 
dM 
dx 



7a 



bj 



l + a"^ M, 
a dM 



|MxHe//| 



1 + a2 M„ dx 



(He// X M)| da; 



wall will move back to its original location. To see this, 
we set 6j = in Eq. (15) and the velocity for H^xt = 
is 



Note that the first term is a pure damping, the second 
and third terms are the result of the spin torque. By 
inserting the effective field Eq. (4), we can express above 
two equations in terms of the wall distortion parameters 
ip{t) and c{t). The expression is algebraically tedious 
and lengthy; we do not write them down here. Instead, 
we would simply consider limiting cases where the spin 
torque is small. A straightforward calculation leads to 
simple but insightful expressions for the rate of the energy 
change and the wall energy, 

dE WfionWoM^ . d^ 

— — = -smificosu}—— (23) 

dt 3 ^ ^ dt ^ ' 

and 

AE = E{bj) - E{Q) = ^L^bj (24) 

where we have defined the inductance of the domain wall 
Lw = 7r/(3/uo7^Wo). These expressions illustrate the en- 
ergy process during the domain wall motion. The total 
spin current induced wall energy is proportional to the 
square of the current density — an analogy to the in- 
ductance of a circuit. 

The stored domain wall energy in the form of domain 
wall distortion can be released once the spin current is 
turned off - an analogy to the electric discharge of in- 
ductance. In fact, we can easily show that the domain 



=(0 



"/AttMs sin ip cos (p 
{l + a^)c{t) 



(25) 



Noted that the initial condition is now ip = p{oo) and 
the initial velocity is +bj/{l + a^). By dropping the first 
two terms in the right hand side of Eq. (14), one can 
integrate (p{t) out from Eq. (14) and it is easy to check 
that 



(22) 



/■oo /"O 
= / Vrecdt = 

Jo J(p( 



dip 



Xn 



(26) 



The domain wall moves to its original position. Again, 
this phenomenon is analogous to the discharge of an in- 
ductance in a usual electronic circuit. We will show in 
the next Section how this feature differs from the domain 
wall motion when an external field is turned off. 



E. Beyond Walker's limit 

Up till now, our analysis is based on Walker's trial 
function. As in the case of field-driving domain wall mo- 
tion, it is essential to assume that the deviation from the 
original Neel wall is small. This places the validity of our 
analytical result to a small current density. In this sub- 
section, we show that when the ciirrcnt density is larger 
than a critical current. Walker's trial function fails to de- 
scribe the domain wall motion. In fact, the domain wall 
will not stop at a sufficient large current and the distor- 
tion of the domain wall can not be simply described by 
two parameters {ip{t) and c{t)) anymore. Fortunately, for 
the experimental interesting, the applied current density 
is usually within the applicability of the Walker's theory. 

To find the critical current density, be, we consider a 
time independent solution by setting ^ = in Eq. (14). 
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By using Eq. (12), we find 

An J Ms sin (p cos ip 



2A 



Maximizing hj with respect to the angle Lp, we have 



' •P'c 



sin cos ipc 



2Hk + 47rMs 
2Hk + 47rMs 



(27) 

(28) 
(29) 



Inserting Eqs. (28) and (29) back into Eq. (27), we find 
the maximum spin torque for a stationary solution is 



when An Ms > Hk, 



AttjMs^ 



(30) 



(31) 



where ^ = ^/AjjAirMf) is the exchange length. Thus 
when bj > be, the stationary solution of Walker fails. 
On the other hand, if the spin torque is smaller than 
be, the maximum distortion shown in Eq. (28) is small 
and thus the Walker's trial function is valid. For a 
Co nanowire, An Ms = 1.8 x 10'' Oe, Hk = 500 Oe, 
7 = 1.9 X 10^(Oe)-is-i, Ms = 14.46 x lO^A/m and 
A = 2.0 X 10~^^J/m, the critical current is about be = 
922m/s or = 2.2 x lO^^A/cm^. This current density 
is at least one to two orders of magnitude larger than 
typical experimental values and thus the trial function 
and its solutions are applicable to most of the experi- 
ments. In Fig. 3, we show the maximum deformation 
of the domain wall; even for the current density as high 
as 2 X 10^'^ A/ cm?, the maximum angle ip is about 8 de- 
grees and the wall only shrinks by 22%. Thus we confirm 
that the trial functions are the good approximation for 
domain wall dynamics at least in one-dimension model 
presented in this section. 



IV. NUMERICAL RESULTS 

Up till now, our description of domain wall dynamics is 
based on the assumption that the out-of-plane angle ip is 
time-dependent but spatially independent. In a typical 
nanowire, the magnetization varies along the direction 
of the wire width, i.e., one may need to go beyond 1-d 
model. Furthermore, a realistic wire is expected to con- 
tain various pinning centers and it would be interesting to 
see the effect of pinning on domain wall motion. In these 
cases, a micromagnetic calculation is required to quanti- 
tatively investigate the domain wall dynamics. Here we 
take an example of a Co nanowire. To make our discus- 
sions concrete, we specify the wire dimension: the thick- 
ness is d = 5 nm, the width is D = 200 nm and the length 
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FIG. 4: The snap shots of the moving magnetic wall at t = 
0.06ns in a Co nanowire. The damping constant a=0.1 and 
the external field He^t = OOe. 



is 1.2 /im. Clearly, this choice preserves one-dimensional 
feature of the wire because D is larger than the domain 
wall width of Co. The reason for choosing this example 
is to make this numerical calculation directly comparable 
with the analytical results obtained in the last section. 
The magnetic wire is divided into rectangular cells whose 
dimension is 2 x 2 x 5nm^. For our film with the large 
width D, the magneto-dipolar interactions between the 
cells can be discarded. The wall dynamics was calculated 
using a fourth order predictor-corrector method with an 
0.8 ps time step. The stationary Neel wall is centered at 
x = at t = 0, a spin current is turned on along x direc- 
tion, and the external field is applied along x direction 
at i > 0. 

In presence of the spin torque, the magnetization of 
domain wall is no more confined in the plane of the layer 
and the original Neel wall is distorted during the motion 
of the domain wall. In Fig. 4, we show snap shots of 
magnetization patterns at i = 0.06 nanoseconds after a 
current is turned at t = 0. For a small current density, 
the shape of the Neel wall is essentially unchanged except 
with a small out of plane magnetization, in consistence 
with our analytical model in previous Section. When 
the applied current is increased to values larger than the 
critical current, the shape of the Neel wall is completely 
destroyed. The magnetization inside the wall has devel- 
oped a significant out-of-plane component, see Fig. 4b, 
and the wall has split into complicated multiple walls. In 
our calculation, the critical spin torque be ~ 1240to/s, 
which approximately agrees with our analytical result 
from Eq. (30) . We now should confine our discussion be- 
low to the case where the current density is smaller than 
the critical current. For a perfect wire, the Neel wall 
moves with a velocity « = immediately after the ap- 
plication of the current. In Fig. 5, we show wall positions 
at different times as well as wall velocity and wall distor- 
tion when wall moves to different positions. The domain 
wall remains to be a quasi-Neel wall and eventually the 
wall completely stops. At a; = 0, the wall velocity is at 
the maximum (vmax = 750to/s). When the wall stops, 
the maximum out-of-plane component of the magnetiza- 
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FIG. 5: (a) The wall positions at several time t = 0, 0.5ns, 
1.0ns, 1.5ns, 2.0ns, 2.5ns and n=3.0ns; (b) the wall velocity 
as a function of the wall position; (c) the maximum out-of- 
plane component of the magnetization as a function of the 
wall position. The damping constant is a = 0.008, the spin 
torque is hj = —750m/s and the external field Hext = OOe. 
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FIG. 6: Domain wall velocity as a function of time, (a) for 
different spin torques at a fixed magnetic field (10 Oe), and (b) 
for different magnetic fields at a fixed spin torque (6j — —500 
m/s). The solid curves represent the analytical results and 
scattered points are the numerical solutions. 



tion (A'/^-componcnt) approaches (\Mz\/Ms = 0.1) and 
the total displacement is Xmax = 312nm. 

It is interesting to take a look at the domain wall dy- 
namics with simultaneous application of the current and 
the external field. As wc already pointed out that the 
initial domain wall velocity is determined by the cur- 
rent but the final terminated velocity is controlled by 
the external field, sec Fig. 6. We confirm here that the 
spin torque alone is unable to move the domain wall 
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FIG. 7: Displacement of the domain wall for several spin 
torques without (a) and with (b) the external field. The solid 
curves represent the solutions of the trial function and scat- 
tered points are the numerical results. 
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FIG. 8: Energy power (a) and net energy gain (b) of the 
domain wall as a function of time. 



over a long distance, see fig. 7(a). This feature is just 

opposite to domain wall motion driven by an external 
field: at an initial application of the magnetic field, the 
wall velocity is small (in fact one can easily show it is 
aHextWo/{l + a^) from Eq. (15)), but it becomes faster 
and faster until it reaches a saturated speed from Eq. (16) 
at Vs = "fHextW {00) /a. The origin of this difference had 
been discussed in great length in the last section and 
we again confirm that these conclusions remain valid for 
the realistic wire. The fact, that the initial domain wall 
velocity is determined by the current but the final termi- 
nated velocity is controlled by the external field, can be 
used to design the high/low-spccd domain wall propaga- 
tion by properly choosing the desired spin current density 
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and the magnetic field. For example, if we are to move 
the domain wall by 75 nanometers, the magnetic field (10 
Oe) alone will take 0.43 nanoseconds. With help of the 
spin torque bj = — 800m/s, the same displacement can 
be made in 0.11 nanoseconds and it takes 0.65 nanosec- 
onds when bj = 300to/s. These featm^es are clearly il- 
lustrated in Fig. 7(b). 

We show the rate of energy input of the domain wall 
during the application of the spin torque in Fig. 8. The 
rate of energy pumping from the spin torque increases at 
initial application of the spin current, then it develops 
a maximum rate. After that, the rate decreases due to 
increase of the damping. Finally, the energy pumping is 
exactly compensated by the damping and the wall motion 
stops. 

It is interesting to compare our calculations with ex- 
periment ^1,] where a domain wall speed about 3m /s for 
a current density 1.2 x IQ^A/crn? was observed. As we 
have concluded that the spin torque alone is unable to 
move the wall over long distance (l//m), the observed 
domain wall velocity over large distance must be from 
two possible scenarios: 1) An additional magnetic field 
generated by the current (Oersted fields) was present in 
the experiment. One just needs a very small magnetic 
field to generate this low velocity domain wall motion. 
In fact, one to two Oersted was sufficient and it was hard 
to rule out such small field in experiment. 2) There is 
a possibility of other spin-current induced torque which 
has different form as Eq. (2). Indeed, if one relaxes the 
adiabatic process that assumes the direction of the spin 
polarization of the current adiabatically follows the mag- 
netization vector in the domain wall, one would gener- 
ate other spin torque. The detail calculation of the spin 
torque from non-adiabatic process is beyond scope of the 
present work; we will defer to a future publication. Other 
experiments [13, (3^, reported that a magnetic field is 
always required to move the wall over a long distance. 
These experiments agree with our general conclusions. 

We end the section by looking at two more interesting 
examples. The first one is that the spin current is ap- 
plied in the form of a pulse. If we applied a rectangular 
pulsed spin current along x direction, what dynamic be- 
haviors do the domain walls have? In Fig. 9, we show the 
evolution of velocity v{t) and displacement x{t) with the 
pulsed spin toque and we compared them with the pulsed 
magnetic field. Suppose the spin torque with the magni- 
tude bj = — 600r7i/s is turned on at t=0 and turned off 
at T = 0.5ns. In the absence of the external field, the 
domain wall moves back to its original location after the 
current is turned off. This is again consistent with our 
analytical result. However, for the field-driving domain 
wall motion, the domain wall continues to move in the 
forward direction after the field is turned off. 

Another example is to model the effect of pinning cen- 
ters in the domain wall motion. In realistic nanowires, 
domain walls are not completely free to move. There are 
various pinning sources such as defects and roughness. 
We introduced a pinning field, i?pi„ = VqxH{C, — |a;|)/C, 




TtTB{ns) Tlme(ns) 



FIG. 9: Velocity v and displacement x for a pulsed spin cur- 
rent (left panels) and for a pulsed magnetic field (right pan- 
els). The damping constant is a = 0.02. 
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bj (m/s) 

FIG. 10: The phase boundary of pinning and depinning do- 
main wall for two types of pinning centers, see text for the 
definition of the pinning parameters Vb and The external 
field is applied along the —x direction and the spin current 
is along -l-a; direction at t = 0. The damping constant is 
a = 0.02. 



where H{y) is the Heaviside step function, ^ > is the 
distance from the domain wall center, x is the position. 
In Fig. 10, we present the result for the critical mag- 
netic field required to move the wall in the presence of 
the defect. The role of the current is to first move the 
domain wall out of the pinning center so that the critical 
magnetic field required to move the wall is smaller. 
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V. CONCLUSIONS 

We have demonstrated in this paper that the domain 
wall motion driven by a spin current has many unique 
features that do not exist in the conventional domain 
wall motion driven by a magnetic field. By solving the 
Landau-Lifshitz-Gilbert equation along with the current 
induced spin torque, we summarize below our main find- 
ings that arc supported by our analytical and numerical 
study: 1) the spin torque alone can move the domain 
wall, but the total displacement is limited, thus the spin 
torque is not capable to move domain walls over a large 
distance; 2) since the speed of the wall is large at an ini- 
tial application of the spin torque, there is an advantage 
of the spin torque over the magnetic field for wall move- 
ment at short distance; 3) the domain wall is capable to 
store current-induced energy and thus the energy process 
of the domain wall is very similar to the charge/discharge 



of a inductance in an electric circuit; 4) the spin torque 
can help to reduce the critical field needed to move do- 
main wall if domain wall pinning centers are present; 5) 
it may be important to extend a spin torque other than 
adiabatic spin torque used here in order to study the con- 
tribution from the non-adiabatic spin torque whose form 
has not yet been investigated. 

In conclusion, a spin torque with the form of Eq. (2) 
is proposed to study the dynamics of the domain wall. 
Due to different roles played by the spin torque and by 
the external field, the current driven domain wall motion 
represent new types of torques. This new spin torque 
creates opportunity for a fast domain wall motion by us- 
ing the combined field and spin torque. The research was 
supported by NSF (ECS-0223568 and DMR-0314456). 

Note added in proof, a paper was published, [? ]atara 
in which the critical current of domain wall motion was 
briefly discussed. 
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